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ABSTRACT 

We investigate the correlation between extinction and H i and CO emission at intermediate 
and high Galactic latitudes (|/7| > 10°) within the footprint of the Xuyi Schmidt Telescope 
Photometric Survey of the Galactic anticentre (XSTPS-GAC) on small and large scales. In 
Paper I (Chen et al. 2014), we present a three-dimensional dust extinction map within the 
footprint of XSTPS-GAC, covering a sky area of over 6,000deg^ at a spatial angular res¬ 
olution of 6 arcmin. In the current work, the map is combined with data from gas tracers, 
including H i data from the Galactic Arecibo L-band Feed Array H i survey and CO data 
from the Planck mission, to constrain the values of dust-to-gas ratio DGR = Ay/N{V{) 
and CO-to-H 2 conversion factor Xco = A^(H 2 )/ITco for the entire GAC footprint exclud¬ 
ing the Galactic plane, as well as for selected star-forming regions (such as the Orion, 
Taurus and Perseus clouds) and a region of diffuse gas in the northern Galactic hemi¬ 
sphere. For the whole GAC footprint, we find DGR = (4.15 ± 0.01) x 10“^^magcm^ and 
^co = (1-72 + 0.03) X 10^°cm“^ (Kkms“')“'. We have also investigated the distribution 
of “CO-dark” gas (DG) within the footprint of GAC and found a linear correlation between 
the DG column density and the V-band extinction: A(DG) ^ 2.2 x 10^* (Ap — A°)cm“^. 
The mass fraction of DG is found to be /dg ~ 0.55 toward the Galactic anticentre, which is 
respectively about 23 and 124 per cent of the atomic and CO-traced molecular gas in the same 
region. This result is consistent with the theoretical work of Papadopoulos et al. but much 
larger than that expected in the H 2 cloud models by Wolfire et al. 
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1 INTRODUCTION 

The interstellar medium (ISM) makes up between 10 to 15 per cent 
of baryonic mass of the Milky Way. About 99 per cent of the mate¬ 
rial is in gas form and the rest is in solid dust grains. Observation of 
the interstellar dust and gas is a primary tools to trace the structure 
and distribution of molecule clouds, where the star formation takes 
place. The dominant molecular gas in the ISM is H 2 , a homonuclear 
diatomic molecule with no permanent dipole moment and allowed 
dipole rotational transitions. The most easily excited H 2 transition 
is the quadrupole S (0): A —> 7/ = 2 —» 0, which has an exci¬ 
tation energy tsEiolk ~ 510 K, more than an order of magnitude 
higher than the typical kinetic temperature ~ 10-50 K of the 
ISM. As a consequence, H 2 is hard to excite and observe, and one 
has to trace the molecular gas and determine the H 2 column den¬ 
sity, A(H 2 ), by other indirect methods. Two of the most commonly 
used tracers are CO gas emission and dust extinction. 
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Heteronuclear diatomic molecule CO is the most abundant 
molecule in the ISM after H 2 . The *^CO 7=1^0 tran¬ 
sition has a REio/k about 5.5 K and can be easily excited even 
in cold molecular clouds and observed from the ground. As a 
consequence, CO has become the primary tracer of the inter¬ 
stellar molecular gas. The molecular abunda nce ratio CO/H 2 is 
about 10“"'. The CO (7 =1—>0) survey of IPame et alj j200lh 
is the most complete survey which covers the Milky Way of 
Galactic latitudes \b\ C 30°. There are a large number of other 
CO surveys covering small sky regions around the individual 
molecular clouds, suc h as the Orion and the Monoceros cloud s 
jHartmann et alJll998l: iMagnani et al] I 2 OOOI : IWilson et al.l l2005h . 
IPlanck Collaboration et al.72013h have extracted the full-sky CO 
emission in three lines, the 7 =1—+0, 7 =2—>1 and 7 =3^2 lines 
at 115, 230, and 345 GHz from the Planck data using a compo¬ 
nent separation method. The *^CO(7 =1—>0) emission at wave¬ 
length A =2.6 mm is then utilized to determine the total col¬ 
umn density of molecular gas assuming a CO-to-H 2 conversion 
factor, the so-called X factor, which relates the surface bright¬ 
ness of CO emission and the column density of molecular gas. 
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^co = ^(H 2 )/Wco (cf. iBolatto et al.|[2013l for a review), where 
Wco is the velocity-integrated brightness temperature of '^CO 
7=1^0. The Xco factor works best for cloud ensemble 
averages. The minimum averaging scale would b e that of a sin¬ 
gle G i ant Molecular Cloud (GMC), --10-100 pc l Dickrnan et al 


200C: [^padoDoulos et alj 


198f: lYoung & Scovillel 1991; iBrvant & Scoville l 19961 ; iRegar 


20020 . The most quoted value of X fac¬ 
tor for the Mi lky Way is Xco ~ 2 x 10^°cm~^ (K kms~‘)"‘ 
jFrerking et 19821 ; iDame et ^l200ll ; iLombardi et al.ll200^ . 

The visual interstellar extinction Ay has been found to have 
a linear relationship with the total column density of hydrogen 
nuclei M(H). The relationship is a cornerstone of the modem as¬ 
trophysics. The classic value of the constant of proportionality, 
i.e. the dust-to-gas ratio, DGR = Av/N(li), derived from the 
measurements of the ultraviolet (UV) absorption of the H i Lycr 
line and the H 2 Lyman band towa rd a small number of early-type 
stars with the Copernicus satellite iSavage et al.ll977l;lBohlin et al 


Il978l). is D GR = 5.34 x 10^22 magcm^. More recently. Iki 


iszl 


( l2014al) and iLisztl j2014bh trace A^(H i) and E{B — V) across the 
sky of Galactic latitudes \b\ = 20 - 60°. They consider data only 
of relatively low extinction, 0.015 < E{B — V) < 0.075 mag, 
for which the hydrogen should mostly be in the form of neutral 
atoms and thus the corrections of A^(H i) for saturation and H 2 
formation are unlikely important. From the data, they obtain a 
DGR value (3.73 x 10~^^ magcm^) that is smaller compared to the 
reference value. The column density of atomic hydrogen, A^(FI i) 
can be determined by measuring the H i 21 cm emission. Thus if 
one can accurately determine the extinction and value of DGR, 
then the H 2 column density can be derived using the relation, 
M(Fl 2 ) = AyjDGR — A^(FIi). There are several successful large 
FI I 21 cm emission surveys of the Gala xy, including the Le iden- 
Argentine-Bonn (LAB) full-sky survey dKalberla et al J 120051) . the 
international Galactic Plane Surv ey (IGPS.lTavlor et alj2003l). the 
Galactic All Sky Survey (GASS. lMcClure-Griffiths et al.l 2003l in 
the southern hemisphere, and the Galactic Are cibo L-band Feed 
Array FI i survey (GALFA-FH. IPeek et al.ll201 lb targeting both the 
Galactic plane and high latitude regions. 

Dust grains can be traced by extin ction in the visible and in 
the near-infrared (IR). IChen et all ( l2014h estimate values of extinc¬ 
tion and distance to over 13 million stars catalogued by the Xuyi 
1.04/1.20 m Schmidt Telesco pe photometr i c survey of t he Galac¬ 
tic an ticentre (XSTPS-GAC; IZhang et all |2013|, I20l4 iLiu et al.l 
I2OI4I) . They derive the extinction by fitting the spectral energy 
distribution (SED) for a sample of 30 million stars that have sup- 
plementar y IR photometry from the Two Micron All Sky Survey 
(2MASS; fskrutskie et al.l 19^ an d the Wide-Field Infrared Sur¬ 
vey Explorer (WISE, Wright et^l2010h . Together with the photo¬ 
metric distance s deduced from th e dereddened XSTPS-GAC opti¬ 
cal photometry. IChen et alj | |2014ll construct a 3D r-band extinction 
map within the GAG footprint at a spatial resolution, depending 
on the latitude, th at varies between 3-9 arcmin. The 3D map of 
IChen et alj ( 1201 4h has also been used to construct a 2D extinction 
map, integrated to a distance of 4kpc, approximately the maximum 
distance that dwarf stars in the sample can trace with sufficiently 
high photometric accuracy. The dust extinction within this distance 
is well constrained. The map has an angular resolution of 6 arcmin 
and is also publicly availablfl In the current work, this 2D extinc¬ 
tion map is used to investigate the dust-to-gas ratio and the X factor 
within the footprint of XSTPS-GAC. 


* http://162.105.156.249/site/ Photometric-Extinctions-and-Distances. 


There is a number of well known large star-forming regions 
within the footprint of XSTPS-GAC, such as the Perseus, Tau¬ 
rus and Orion clouds. They form parts of the Gould Belt. Those 
clouds are in the souther n Galactic plane, l ocated at relatively small 
distances {d < 800pc; Ichen et alj|2014h . We study the correla¬ 
tion between the interstellar dust and gas for these regions. Within 
these regions there exist diffuse, non self-gravitating clouds as 
well as dense, massive, self-gravitating, star-forming clouds like 
the Orion. Several studies on the dust and gas correlat i on hav e 
been carried out for tho se in dividual cloud s. Pineda et alJ l l2008ah . 
[Goldsmith et alj l l2008h and Pineda^eLalJ ijoiol ) study th e Tau¬ 
rus molecular cloud, while [Pineda et alj ( l2Q08bll ! iLee et all ilOl'j) 
and Lee e t alj ll20 14h focus on t he Pe rseus cloud a nd Digel et ^ 
( I 1999 I 1 . Ackermann et alj ( 1201 2h and iRipple et alj ( 1201 3n on the 


Orion cloud. 

One can measure the amount and distribution of the atomic 
and molecular gas as well as the dust grains by combining var¬ 
ious tracers. In this paper, we combine the dust extinction to- 
geth er with gas distrib ution traced by H i data from the GALFA- 
HI jPeek et al. 2011 ) and CO data from the Planck mission 
dPlanck Collaboration et aljl2013h . to constrain the values of DGR 
and Xco at a high angular resolution of 6 arcmin. The exis¬ 
tence of molecular gas of low C O abundance has been noted 
previously in theoretica l models dvan Dishoeck & BlackI 1 19881 : 
IPapadopoulos et a l]|2002|). in obser vations of diffuse gas and high- 
latitude c louds ( Lada&Blitz 19881) and in observations of irregular 
galaxies dMadden et alJll997 ). The values of DGR and Xco deter¬ 
mined in the curren t work allow us to further trace this so-called 
“CO-dar k gas” (DG; Reach et all 1994 ; Meverdierks & HeithauserJ 


ll996l:lGrenier et al.l200.5l:lAbdo et al.l201(lh at both small and large 

scales. We are not the first one on this topic of course. For example, 
[Planck Collaboration et al.l ( 1201 ih has obtained an all-sky map of 
dust optical depth and compared it with the observed distribution 
of gas column density. The results are used to constrain the correla¬ 
tion between the dust and gas, {jo/NyiY‘^, as well as the values of 
Yco> and to construct a distribution map of DG at high an d interme¬ 
diate Galactic latitudes (|Zi| > 10°). |Paradis et alj ( l2012h carry out 
a similar study using the extinction data deduced from the colour 
excess of the 2MASS photometry instea d of those derived from 
the far-IR optical depth. Both the work of IParadis et al.l (1201 2h and 
IPlanck Collaboration et al.l ( 1201 ih are limite d by the angular reso - 
lution of the LAB 21 cm H i emission survey dKalberla et al.l20()l) . 
which is 36 arcmin. In the current work we have adopted data from 
the GALFA-HI survey. As a consequence, the angular resolution 
has been improved significantly. 

The paper is structured as the following. In Section 2 we 
present the 2D extinction map integrated to a distance of 4kpc. 
Section 3 describes the data of gas tracers used in the current anal¬ 
ysis. The correlations between extinction and gas content are in¬ 
vestigated in Section 4 with the main results discussed in Section 5. 
We summarize in Section 6. 


2 THE EXTINCTION MAP 

The p hotometric extinction and distance catalogue of IChen et al.l 
d2014l) is first used to construct a 2D extinction map integrated to a 
distance of 4kpc for the footprint of XSTPS-GAC. The catalogue 
contains more than 13 million stars with estimates of r-band ex¬ 
tinction and photometric distance. To ensure good photometric ac¬ 
curacy, only objects flagged with a quality flag ‘A’, i.e. with average 
photometric errors less than 0.05 mag for all the eight bands used 
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to derive the extinction and distance (g, r, and i from the XSTPS- 
GAC, J, H and from the 2M ASS, and Wl and W2 from the 
WISE), are included ISamnle A of ichen et al.ll2oI3l . Possible can¬ 
didates of giants, flagged as ‘G’, are excluded because their dis¬ 
tances may have been grossly underestimated. This leaves just over 
7 million stars and they are divided into sub-fields of size 6 x 
barcmin^. For each sub-field, a sliding window of depth 450pc, 
stepping by steps of 150 pc, is used to obtain the median value of 
extinction for each distance bin. The extinction as a function of 
distance, A'fd) = f{d), where i is the index of sub-fields, is then 
obtained after applying a boxcar smoothing of width 450 pc. The 
extinction versus distance relations of all sub-fields (sightlines) are 
used to construct a 3D extinction map, which is then integrated out 
to a distance ofd = 4 kpc to yield a 2D extinction map to be used in 
the current analysis. Most of the sample stars are within a distance 
of 4 kpc and the extinction out to this distance is well constrained. 
The 4 kpc integrated map can be treated as a lower limit of extinc¬ 
tion map along the sightlines. 

Fig. U presents the 4 kpc dust extinction map thus derived. 
It has a spatial resolution of barcmin. The resolution is chosen 
to ensure that most of the subfields have a sufficient number of 
stars (> 10) to derive a robust relation of extinction as a func¬ 
tion of distance. For a small fraction of sub-fields (less than 1.5 
per cent) of high Galactic latitudes or suffering from high extinc¬ 
tion, there are not enough numbers of stars of high photometry ac¬ 
curacies available. For those fields, we have adopted some stars 
of lo wer accuracy from the catalogue (Sample C of IChen et ^ 
l2014h . Those sub-fields are flagged by ‘C’. The V-b and extinc¬ 
tion A y is converted from using the extinction law of lYuan et al.l 
( l2013l) , A„ = 1.172A,, assuming a total-to-selective extinction ra¬ 
tio Ry = 3.1. The uncertainties of the map are estimated from 
the stars located within the distance bin of width 450 pc centered 
at 4 kpc and in a width of 450 pc. The typical uncertainties are 
O.lSmag in Ay, or about 0.06mag in E{B — V). The uncertain¬ 
ties are mainly determined by the stellar density and depth of the 
sample. Under favorable conditions, the uncertainties could be as 
low as 0.04 mag in Ay. Some “stripe-like” artifacts can be seen on 
the map. They are caused by the nightly varying survey depth and 
photometric accuracy of the XSTPS-GAC. A notable empty patch 
around I ~160° and b ~20° is due to the abnormally large photo¬ 
metric errors of XSTPS-GAC in that area. 

Chen et al. (2014) have com pared their Sample A me- 
dian extinction map with those of ISchlegel et alJ i ll9981) and 
iFroebrich et al.l | |2007|) . and found some systematic offsets amongst 
them. Parts of the differences are clearly caused by the difference 
in distance traced by the three maps [2-3 kpc for Sample A me¬ 
dian extinction map of Chen et ah , essentially infi n ite for the SFD 
map and about 1 kpc for that of IFroebrich et alj | |2007|) . respec¬ 
tively]. The 4 kpc extinction map presented here should have in¬ 
cluded almost all the dust extinction alone the sightlines for sub¬ 
fields of intermediate and high Galactic latitudes {\b\ > 10°), thus 
should be directly comparable with that of SFD. In addition, it 
should also be comparable wit h the 4.5 kpc redde ning map newly 
derived from t he Pan STARRS jKaiser et alj|2010l) photometry by 
ISchlaflv et al.l j2014|) . Fig. [2] compares values of e xtinction from 
our 4 kpc extinction map with those of SFD and ISchlaflv et al.l 
( l2014h for sightlines of Galactic latitude \b\ > 10°. Values of 
Ay are converted to E[B — V) assuming Ry = 3.1. Except for 
some outliers, our v alues agrees well with both those of SFD and 
ISchlaflv et al.l |2014|) for E{B — V) < 0.8 mag. In deriving th e ex¬ 
tinction, both our analyses and those of l^hlaflv et alj (l2014h have 
made the assumption that all reddening values should be no less 
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Figure 1. A 2D extinction map, integrated to a distance of 4 kpc for the 
footprint of XSTPS-GAC. The map has an angular resolution of 6 arcmin. 


than zero. The assumption could introduce some biases, in partic¬ 
ular for sightlines of very low, close to zero, extinction. I n princi¬ 
ple, both the best SED fit al gorithm of Chen et al. ll l20Ti) as well 
as the Bayesian approach of iGreen et al.l J 2 OI 4 I) could lead to sys¬ 
tematically oyeHestimatedvalues_of_extinctonJb£jcgions_of2ow 
extinctio n dSerryet al.l2012LIChen et alj2014l : [Schlaflv et al.l2014l : 

I Green et~^ .1 120141) . Photometric errors obviously contribute to the 

uncertainties. Other potential sources of error include the colour¬ 
reddening degeneracy. Those systematics might be the cause of the 
significant non-Gaussianity seen in the distribution of differences of 
values as given by our map and those of SFD (Fig. 2). By contrast, 
the distribution of di fferences of values between ours and those of 
ISchlaflv et al.l (1201 4b follows closely a Gaussian one, reflecting the 
fact that both use a similar method to estimate the extinction. 

For high extinction regions [E{B —V)> 0.8 mag], SFD red- 
dening values are sy stematically higher than ours while those of 
ISchlaflv etai]|20l4 map are lower. Most of those highly reddened 
regions for which our estimates of reddening are significant smaller 
than given by the SFD map are actually from low Galactic lati¬ 
tudes. In fact, if we restrict the plots of Fig. 2 to absolute latitudes 
higher than 20° (\b\ > 20°), then the systematic differences be¬ 
tween ours and the SFD values disappear. Thus it is very likely, 
for those highly reddened, low latitude regions, the integration dis¬ 
tance of 4 kpc used to derive our current 2D extinction map is not 
deep enough. The s ystematic lower values given by the map of 
ISchlaflv et al.l | |2014|) compared to ours for those highly extincted 
regions is much more difficult to understand, as both analyses use 
similar methods, although the results are based on different data 
set. The integration of the Schlafly et al. map is 500 pc further than 
ours. But instead of giving higher values, the map of Schlafly et 
al. act ually yields lower values compared to ours. ISchlaflv et al.l 
(l2014h notice that when their map is compared to that of SFD, their 
map yields systematic lower values even for regi ons of latitudes 
\b\ > 30°. It is possible that lSchlaflv et alj ( |2014|) may have sys¬ 
tematically underestimated the reddening values for some of the 
highly extincted regions. For the whole XSTPS-GAC footprint of 
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highest angular resolutions and best signal-to-noise ratios (SNRs) 
are employed in this work. 
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3.1 HI data 

We use the H i data from the GALFA-HI survey to trace the content 
of atomic gas. The GALFA-HI survey uses Arecibo L-band Feed 
Array, a seven-beam array of receivers mounted at the focal plane 
of the 305 m Arecibo telescope , to map the H i 21 cm emission 
in the Galaxy dPeek et al.ll201lll . The survey aims to cover both 
the Galactic plane as well as the high Galactic latitudes with an 
angular resolution of 3.4 arc min. The GALFA -HI first data release 
(DRl) has been presented bv IPeek et alJdTOllh . The DRl data cov¬ 
ers 7520 deg^ of the sky, produced from 3046 hr worth of data ob¬ 
tained for 12 individual projects. The local standard of rest (LSR) 
velocity of the emission ranges from —700 to 3-700kms“'. The 
rms noise in a lkms“' channel ranges from 140 mK to 60 mK, 
with a median of 80 mK. The DRl data can be downloaded from 
https://purcelLssl.berkeley.edu/. 

The GALFA-HI DRl data does not cover the entire footprint 
of XSTPS-GAC, but only the region of 3h < RA < 9h and 
0 < Dec < 39 deg. This unfortunately restricts the sky cov¬ 
erage of the current work. We use the data of velocity resolution 
0.8kms“'. The H i column density A(H i) is derived from the 
profile of brightness temperature after applying a small correction 
for the optical depth effects assuming a spin temperature of 145 K, 
as adopted bv iLisztl (I2014al l9). The temperature c orresponds to the 
mean ratio of H i 21 cm emission and absorption dLiszt et alj2010ll 
measured by the Millen nium H i Emission-Absorption Survey of 
iHeiles & Trol^ ( l2003l) . Values of A(H i) derived for the XSTPS- 
GAC footprint range from 1.5 xlO^^to 1.2 X 10^^ cm~^, with a me¬ 
dian of 1.4 X 10^* cm^^. The distribution of A(H i) thus deduced is 
shown in the left panel of Fig.|3] 


Figure 2. Comparisons of extinction values as given by our map and those 
given by the SFD map (top two panels), and by the map of Sclafly et al. 
(2014; lower two panels). The red pluses in the first and third panels rep¬ 
resent median values in the individual bins. A black straight line denoting 
complete equality is also overplotted to guide the eyes. In the second and 
fourth panels, the black curve is a Gaussian fit to the distribution of differ¬ 
ences of values. 


\b\ > 10°, our map gives on aver age only marg i nally larger val¬ 
ues compared to those of SFD and lSchlaflv et al.[ ( l2014h . by 0.013 
and 0.015 mag, respectively. The scatters of differences between the 
values are about 0.06 mag in E{B — V), comparable to the typical 
uncertainties estimated for our map (0.18 mag in Ay). 


3 GAS TRACERS 

The interstellar gas can either be atomic, molecular or ionized. 
Those three phases are commonly traced by the H i 21 cm emis¬ 
sion, the molecular CO pure rotational transitions and by the H i 
Ha recombination line, respectively. The ionized gas can also be 
traced by the free-free emission of thermal electrons. The contri¬ 
bution of the ionized gas, which are mostly restricted to the indi¬ 
vidual, isolated H it regions, to the total gas content is neglected in 
the current study. Only most recently published data that have the 


3.2 CO data 


We use the CO-integr at ed in tensity image produced by 
IPlanck Collaboration et al.l ( 1201 3h . which is based on the 
Planck first 15.5 month survey data (two full-sky scans) and 
uses the full-sky maps of the nine Planck frequency bands, 
and also the 100, 217 and 353 GHz full-sky bolometer maps. 
IPlanck Collaboration et al.l ( 1201 3h provide all sky maps of the 
velocity-integrated emission of the CO 7=1^0, 7 = 2—>1 
and 7 = 3 —> 2 transitions. They produce three types of maps. We 
choose the Type 3 map generated from the multi-line approach. 
The approach assumes fixed CO line ratios in order to achieve 
the highest possible signal-to-noise ratios, and has an angular 
resolution of ~ 5.5 arcmin. The standard deviation of the Type 3 
CO(l^O) map, when used at an angular resolution of 15 arcmin, 
is typically 0.16 K, compared to 1.77 and 0.45 K for the Type 1 
and 2 maps, respectively. For comparison, the CO survey of iDamel 
(I 2 OI ih has a typical uncertainty of 0.6 K. We choose the Planck 
Type 3 CO (1-^0) map because of its high angular resolution 
and high SNR. Type 3 map is al so more suitable for intermediate 
and high-Galactic latitudes (cf. [Planck Collaboration et alj 1201 3l 
for more detail). There is some *^CO contamination in the Planck 
CO(l^O) ma p. To correct for this, we have div ided the data by a 
factor of 1.16 dPlanck Collaboration et al]l2013h. 


IPame et^ ( 12001 ) combine data of CO surveys of the Galac¬ 
tic plane and local molecular clouds and produce a composite map 
of an angular resolution of 8.4 arcmin. The data were obtained with 
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Figure 3. Contours of A?(H l) (left panel) and Wco (right panel) in the footprint of XSTPS-GAC. The grey-scale background image illustrates the 4kpc 
extinction map. The squares denote different regions analyzed in the current work. The maps have an angular resolution of 6 arcmin. 


the 1.2 m telescope of the Harv ard-Smithsonian C enter for Astro¬ 
physics. The data published by IPame et^ (1200 ih are mainly for 
the Galactic plane. Data of high Galac tic latitudes {\b\ > 20°) are 
not available. In this paper, the data of lPame et al.l ( 200 ih are only 
used as a reference to check the data we actually use. 

The Planck Type 3 CO-integrated intensity image is presented 
in the right panel of Fig. (3] We only show the footprint where data 
of Ay, A(H I) and Wco are all available. It shows excellent agree¬ 
ment with the dust extinction map. In this paper, we first bring all 
input data sets to a common angular resolution determined by the 
lowest angular resolution data set, which is 6 arcmin of the extinc¬ 
tion map. This is done by performing a Gaussian kernel smooth. 
Thanks to the high resolution, we can analysis the correlation be¬ 
tween dust and gas at intermediate Galactic latitudes of the GAC 
for the individual molecular clouds, such as the Orion, Taurus and 
Perseus cloud. We divide the clouds into eight different regions. 
Their positions and ranges on the sky are marked by the squares 
in Fig. [3] together with a region of diffuse gas without any obvi¬ 
ous molecular clouds in the northern Galactic hemisphere, noted as 
“North” in the current work. 


4 EXTINCTION AND GAS CORRELATION 

The observed visual extinction Ay can be modelled by contribu¬ 
tions from the atomic gas and the CO-traced molecular gas, as 

A7‘' = DGRiNiU i) + 2Xco Wco) + A“ , (1) 


fitting strategy similar to that of IParadis et al.l |2oT2) • The fitting 
process can be described as following. 


(i) We first assume that there is no DG in all subfields (pix¬ 
els). We then select pixels where the CO emission is low (Wco < 
0.2Kkms”*). The column density of molecules of those pixels 
A(H 2 ) must be less than 2.0 x 10*^ cm“^ if one assumes that 
Xco = 1 xl0^®cm ^(Kkms *) *. Given that the minimum col¬ 
umn density of atomic gas within the XSTPS-GAC footprint is 
1.5 X 10“® cm~^ (see Section 3.1), the gas column density of those 
pixels can be considered as being contributed by atomic gas only 
[A(H)=A(H I)]. Under such assumptions, we then derive the best- 
fit values of DGR and Ay of those pixels. 

(ii) The best-fit values of DGR and Ay obtained from the last 
step are then applied to all pixels, and this time the value of 
Xco is fitted as a free parameter. We have smoothed the data 
to an angular resolution of 2°, corresponds to about 140 pc at 
a distance of 4kpc and 15 pc at a distance of 400 pc (where 
most of the molecular clouds analyzed in the current work, the 
Orion, Taurus and Perseus clouds locate within) in this step 
to obtain the statistical cloud-ensemble-average values of Xco 
toickman^et jI]|lj8j;jYoun g & Scovillelfl99ll : IBrvant & ScoviIi3 
1 19961: Papadopoulos et al.ll2002h. 

(hi) The above derived values of DGR, A® and Xco are applied 
to Eq. (1) to calculate values of A™'^ of all pixels. h)f value is then 
calculated = 2, (Ay/ — Ayf), where i is the pixel index. 

(iv) The column density of DG is then calculated from formula. 


where A® is a constant to account for possible measurement offset 
of the extinction map and/or possible contribution to the extinc¬ 
tion from the ionized gas. The total gas column density is given by, 
A(H) = A(H i) -I- 2A(H2). We constrain the three free parameters, 
DGR, Xco and A®, using the visual extinction, H i column density 
and CO integrated line intensity maps described in the last Section. 
Considering the high resolution of the maps and the large uncer¬ 
tainties of extinction for low extinction regions, we have adopted a 


A(DG) = {Ay - A’tffjDGR. (2) 

(v) We then go back to the first step, but this time use the DG 
distribution derived from the previous step. We select pixels that 
both have low CO emission (IVco < 0.2 K km s“*) and no DG con¬ 
tent [A(DG) sS 0cm~^]. The fields are used to update the best-fit 
values of DGR and A® and the process is iterated until a minimum 
value is reached . 
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presented in Fig. [S] Clearly, different clouds seem to have differ¬ 
ent dust and gas properties. The values of DGR deduced for the 
individual regions of molecular clouds range from 3.24 x 10“^^ to 
7.74 X 10~^^ mag cm^. For the diffuse gas region “North”, the value 
is 1.38 X lO'^^magcm^. The values of Xco for different clouds 
range from 0.84 x 10^° to 2.39 x 10^® cm~^ (Kkms“')~*- The 
correlation between Ay and A^(H) is quite good for the individual 
clouds. In the “North” region, the lack of high CO content pixels 
makes it very difficult to fit Xco- The dust extinction in this region 
is also low (most of the pixels have Ay < I mag), leading to a 
value of DGR, much smaller than those found for the regions of 
molecular clouds. 


5 DISCUSSION 

5.1 Comparison with previous work 


Figure 4. Correlation between the visual extinction Ay and the total column 
density of hydrogen, W(H) = N{Hi) + 2N{H2), for the whole XSTPS- 
GAC footprint of |1)| > 10° that excludes the Galactic plane. The blue con¬ 
tours represent the density of pixels on a logarithmic scale. The red pluses 
are median values for the individual A(H) bins. The red solid line is a fit to 
all the data points. 


In the current work, we have omitted the Galactic plane which 
shows very strong atomic gas emission and where the contribution 
of ionized gas can not be ignored. Thus only data from Galactic 
latitudes |i>| > 10° are used throughout the paper. Since we have 
excluded the Galactic-plane clouds, those included are representa¬ 
tive of the diffuse cloud component of the Galaxy. Usually a mini¬ 
mum is reached after the first iteration that assumes there is no 
DG component at all pixels. We list the results in Table 1 for both 
the whole XSTPS-GAC footprint of \b\ > 10° and for the individ¬ 
ual selected regions. Errors of the resultant parameters are derived 
using a bootstrap method. For each parameter, we apply the algo¬ 
rithm to 1,000 randomly selected samples and derive the parame¬ 
ters of those samples. The resultant parameters follow a Gaussian 
distribution. The variance of the Gaussian distribution is taken as 
the error of the corresponding parameter. 

We present in Fig. |4] the correlation between visual extinc¬ 
tion Ay and total column density of hydrogen, Ai(H) = N(H i) -I- 
IXcoWco, for all XSTPS-GAC pixels of \b\ > 10°, excluding the 
Galactic plane. Despite the large scatters, in particular for pixels of 
low extinction, the data exhibit a significant linear correlation be¬ 
tween the two quantities. The data yield a best-fit value of DGR of 
(4.15±0.01) x 10“^^ mag cm^. The best-fit value of Xco is found 
to be (1.72 ±0.01) x 10“ cm'^ (Kkms-‘)-‘. 

The southern Galactic hemisphere covered by the XSTPS- 
GAC survey {b < —10°) include several prominent giant molec¬ 
ular clouds, including the Orion, Perseus, Taurus clouds and the 
Taurus cloud extension regions. We have earmarked eight regions 
and examine the gas and dust correlation. We have also applied our 
algorithm to a region of diffuse gas in the northern Galactic hemi¬ 
sphere, earmarked as “North”. The ranges, mean extinction, mean 
A(H i) and mean Wco of these regions are listed in Table 2. The 
positions of those regions are marked by squares in Fig[3] The CO 
and extinction maps show remarkably similar features, although 
those features are not so obvious in the H i map. For those individ¬ 
ual regions, the resultant parameters are also listed in Table 1 and 
the relevant data showing the correlation between Ay and A(H) are 


The average value of DGR deduced here for the XSTPS-GAC 
footprint that excludes the Galactic plane, (4.15 ± 0.01) x 
10 ~^^magcm^, is in good agreement with the results of pre¬ 
vious work. It is slightly sm aller than th e mos t cited value, 
5.3 X 10““ magcm^, derived bv ISohlin et akl (Il978l) from the UV 
and optical absorption line measurements toward hot stars. It is 
also slightly smaller than values of (6.0 — 6.5) x 10“^^ magcm^ 
derived fro m the reddening and the Lya absorption toward early - 
type stars dShull & vanSteenl^gl 1 19851 : IPinlas & Savagel 119941) . 


On the other hand, Lisztl 1 2014at) finds an even smaller value of 
3.7 X 10““ magcm^ from the SFD map combined with 21 cm H i 
measurements. The values of DGR for the individual cloud regions 
vary by a factor of two, from 3.2 to 7.7 x 10“^^ mag cm^. They are 
also c onsistent with values from the previous work. IParadis et al.1 
( l2012h obtain an average value of 5.46 x 10“^^magcm^ for the 
whole region of the outer Galaxy (|/| > 70°) and values of 
(3.4 — 8.7) X 10““ magcm^ for individual regions, in good agree¬ 
ment with our results. A variety of factors may lead to variations of 
dust-to-gas ratios amongst different regions, including differences 
in dust content and metal abun dance. The dust extinct ion law of dif¬ 
ferent regions may also differ. iFroebrich et alj (l2007h discuss vari¬ 
ations of the extinction law amongst individual clouds in the GAG 
area. In fact, we notice that the value of DGR decreases with in¬ 
creasing power-law index of extinction law in the near-IR. The 
index fi is known to corr elate with the dust properties, including 
the grain size distribution dPrain 3120031) . 

The value of CO to H 2 conversion factor, Xco. deter¬ 
mined hitherto differs from one analysis to another. Our best-fit 
value of 1.72 x 10“ cm“^ (Kkms“*)“' for the whole XSTPS- 
GAC foot print of |fc| > 10° is close to the Galactic aver¬ 

age value dStrong & MattoxI Il996h . The value is slightly lower 
but nevertheless quite close to the value of (1.8 ± 0.03) x 
10 “ c m“^ (Kkms“‘)“', derived for Galactic regions of |fo| > 
5° by IPame et ^ ( 1200 ih based on extinction data from SFD, 
H I data from the Dw i ngeloo-Leiden H i survey and CO data 
from IPame et^ ( 1200 ih . iGlover & Mac Lo^ ( 1201 ill demonstrate 
by numerical simulations that 2fco falls off with extinction as 
Xco QC for clouds of mean extinction smaller than 3mag. 

However, our current results do not show that kind of relation. 
In addition, our derived values of Xco are significantly smaller 
than p r edicted by the num erical simulations of iGlover & Mac Low! 
( I 2 OI ih . IPolk et aP ( Il988h consider a simple two-component cloud 
model which consists of GMC cores and lower column density 
(“diffuse”) molecular gas. The conversion ratios for the dense 
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Table 1. Dust to gas ratio DGR and CO-to-H 2 conversion factor Xco derived for individual regions toward the Galactic anticentre. 


Region 

DGR 

(10“22 magcm^) 

Xco 

(mag) (10^® cm~^ (Kkms“*)“') 

Mh(DG) 

Mh(DG) 

/do 

V(DG) 

(mag) 

Mh(Hi) 

Mh(H“) 

(Av-A°) 

10^^ cm~^ mag~^ 

All of \b\ > 10° 

4.15 + 0.01 

0.06 + 0.002 

1.72+0.01 

0.23 

1.24 

0.55 



Orion 1 

7.74 + 0.21 

-0.75 ± 0.05 

1.54+0.01 

0.09 

0.45 

0.31 

2.15+0.07 

1.10+0.05 

Orion 2 

5.53 + 0.08 

-0.29 + 0.02 

1.34+0.03 

0.10 

1.43 

0.59 

1.69+0.07 

0.97+0.06 

Taurus 

3.65 + 0.10 

0.34 + 0.02 

1.38+0.01 

0.22 

0.31 

0.24 

2.75+0.12 

0.91+0.04 

Taums El 

6.42 + 0.10 

-0.15 + 0.02 

0.84+0.01 

0.11 

1.43 

0.59 

2.28+0.02 

0.76+0.04 

Taums E2 

5.62 + 0.08 

-0.14 + 0.01 

1.69+0.04 

0.11 

4.01 

0.80 

1.89+0.07 

0.64+0.05 

Taurus E3 

6.33 + 0.15 

-0.13 + 0.02 

1.41+0.02 

0.13 

3.42 

0.77 

2.28+0.03 

0.57+0.04 

Perseus 1 

5.01 + 0.08 

0.10 + 0.01 

1.04+0.01 

0.18 

0.89 

0.47 

2.10+0.04 

0.59+0.05 

Perseus 2 

3.24 + 0.28 

0.22 + 0.03 

2.39+0.01 

0.22 

1.05 

0.52 

2.36+0.04 

0.55+0.04 

North 

1.38 + 0.02 

0.20 + 0.002 


0.95 



2.40+0.03 

0.22+0.04 











Table 2. Individual regions studied in the current work. 


Region 

longitude 

latitude 

Ay (mag) 

N{H l) (102°cm-2) 

Wco (Kkms-*) 

Orion 1 

211°< 1 < 201° 

-17° <b< -10° 

1.44 

23.93 

1.61 

Orion 2 

201°< / < 188° 

-18° <b < -10° 

1.04 

22.31 

0.56 

Taurus 

178° < / < 164° 

-19° <b< -10° 

1.38 

16.47 

4.22 

Taurus El 

189°< / < 171° 

-31° <b < -17° 

1.02 

16.78 

0.75 

Taurus E2 

182°< 1 < 165° 

-45° <b< -32° 

0.70 

14.62 

0.12 

Taurus E3 

197°< 1 < 183° 

-41° <b< -28° 

0.66 

11.96 

0.16 

Perseus 1 

165°< 1 < 150° 

-25° < 6 < -10° 

0.93 

13.35 

1.30 

Perseus 2 

171°< 1 < 152° 

-32° <b < -25° 

0.65 

11.49 

0.49 

North 

220°< 1 < 175° 

10°< < 30° 

0.31 

7.70 

0.00 


cloud cores and the diffuse clouds are different, estimated at ~ 
4 X 10^° and 0.5 x 10^° cm^^ (Kkms“*)“*, respectively. Thus 
the measured values of Jlfco in the current work (around ~ 1.5 x 
10 ^° cm~^ (Kk ms~‘)~‘ ) a r e clos er to those of the diffuse clouds 
of the work by IPolk et aP ( Il988h than the dense ones, indicat¬ 
ing that the clouds in our sa mple are dominated by the diffuse 
phases. IPineda et alj ( l2008ah find that Xco varies from 0.9 to 
3.0 X 10^° c m~^ (Kkm s~‘)~‘ among st several regions of the 
Perseus cloud. lAckermann et al.l ( 120120 find that for the Orion 
cloud, Xco ~ (1 .3 ~ 2.3) x 10^® cm~^ (Kkms“‘)“*. Finally, 
lAbdo et al.l ( I 2 OIOI 1 report a significant variations ranging from 
~0.9x to ~ 1.9 X 10^® cm~^ (Kkm s“*)“', from the Gould Belt to 
the Perseus arm. All those findings are all in good agreement with 
our results. 


5.2 Validation of the gas tracers 

Both IPlanck Collaboration et al.l ( 1201 ih and IParadis et al.l ( 1201 2h 
assume that the 21 cm line is optical th i n whe n calculating the col¬ 
umn density of atomic gas. lFukui et al.l(l2014h infers that the DG in 
the Galaxy is dominated by optically thick and cool H i gas, which 
implies that the average density of H i is two times higher than 
that derived assuming optically-thin in the local interstellar space. 
In this work we have applied a first-order opacity correction, as¬ 
suming a constant for the emission. The need to correct for the 
opacity effects in the 21 cm line profiles ultimately limits the accu¬ 
racy of A1(H) determinations. We find that for most directions to¬ 
ward the GAC but out of the Galactic plane, the column densities of 
neutral hydrogen derived are systematically higher than those de¬ 
duced assuming optically thin. A comparison of H i column densi- 



N(HI) T,= 145 K (10^'' cm“^) 


Figure 6. Left panel: Comparison of atomic hydrogen column densities de¬ 
rived assuming a constant spin temperature of T, = 145 K and those de¬ 
duced assuming optically thin emissio n. Right panel: Co mparison of the 
corrected Planck CO data and those of lPame et al.l j200ll) . The red pluses 
are median values for the individual bins and the black solid lines represent 
the line of equal. 


ties derived assuming optically thin and those assuming a constant 
Tj = 145 K is presented in the left panel of Fig.The latter values 
are 1.1-1.3 times higher than their optically thin counterpart. 

In the ri ght panel of Fig. [3 we compare the Planck CO data 
with those of TPame et al.l (120011). The Planck Type 3 CO data are 
in good agreement with those of IPame et al.l ( l200]f) . after correct¬ 
ing for a constant slope of 1.16. We have investigated the parame¬ 
ters fo r the overlapping subfields using the CO data of IPame et al.l 
( I 2 OOII) instead of those Planck. It yields consistent results very sim¬ 
ilar to those presented here. 
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Figure 5. Same as Fig.|4]but for the individual regions. The name of the region is marked in each panel. 


5.3 The CO dark gas 


There is a clear excess of Ay as given by the extinction map 
compared to the simulated A"”'' for regions of Ai(H) > 1.2 x 
10^' cm'2 (Ay > 0.6 mag, see Fig. 0. The excess corresponds to 
the p r esence of CO dark gas, i.e. DG jPlanck Collaboration et al.l 
|201l|;|Paradis et alJ|2012|) . The DG could be interstellar gas that 
exists in the form of F12 along with C i and C ii, but con¬ 
tains little or no CO. The column density of DG, A'(DG), 
can be determined using Eq. (2). The total mass of hydrogen 
of the DG component, M(DG), the mass of atomic hydrogen, 
M(Hi), and the mass of t he CO-traced molecular ga s , 
can be estimated fo llowing [Planck Collaboration et all ( 1201 ih and 
IParadis et al.l ( l2012h : M(DG/H i/H“)cx ^ A'(DG/H i/H“). We 
have assumed that all gas component have the same spatial extent. 
We have M(DG)/M(Hi) = 0.23 and M(DG)/M(H“) = 1.24, 
respectively. The values of those quantities thus derived are pre¬ 
sented in Table 1. Our estimates of DG ma ss are lower than those 
found b vIPlanck Collaboration et ^ ( 1201 ih but are consistent with 
those of IParadis et al.72012h . The estimated DG mass is about 19 


per ce nt of the total gas mass, ve r y clos e to the value of 22 per 

and the value of 16 per 


Planck Collaboration et alj 


•ycJose 

OOllh 


Paradis et al. I (l2012h . The small differences amongst those 


cent of 
cent of 

estimates are however probably insignificant owing to the large un¬ 
certainties involved in the analyses. On the other hand, the “close” 
agreement between the results from different analyses, which uti¬ 
lize different data set of differing spatial resolutions and sky cover¬ 
age probably indicate a “true universal” fraction of DG mass rela¬ 
tive to that of all gas, ~ 20 per cent on large Galactic scale. 


The spatial distribution of DG derived from the current work 
is shown in the upper panel of Fig. |3 for the entire XSTPS-GAC 
footprint where the GALFA-FII data are available. Regions where 
VFco > lKkms ~‘ have been excluded. The distrib ution agrees 
well with that of [Planck Collaboration et al. I j2onh and clearly 
shows that the DG is distributed mainly around molecular clouds, 
such as the Taurus, Perseus and Orion clouds and exhibits features 
similar to the dust distribution. Significant amount of DG is found 
south of the Galactic plane within the footprint of XSTPS-GAC 
but little in the north. If we compare the mass of DG to those of H i 
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and CO amongst different clouds, we see significant variations. For 
example, the ratio of M(DG)/Af(H i) varies from a value of 0.09 
for Orion 1 to 0.22 for the Taurus cloud and to 0.95 for the North 
region. Similarly, the ratio of M{DG)/M{iif') increases from a 
value of 0.31 for the Taurus cloud to 1.43 for Taurus El/Orion 
2 and to 3.42 and 4.01 respectively for Taurus E3 and E2. The 
value of M(DG)/M(H 2 °) for the North region is not listed in Ta¬ 
ble 1 given the extremely large uncertainty as a consequence of 
the lack of detection of CO emission in this region. All the ra¬ 
tios quoted above are however highly dependent on the values of 
DGR and Xco derived. The mass fraction of DG, /dg, is defined 
as /dg = M(DG)/[M(DG) + M(H“)]. Its values for the individ¬ 
ual regions are also given in Table 1. For the entire XSTPS-GAC 
region of |Z)| > 10°, we fin d /dg = 0.55, in good agreement with 
the value of 0.54 found by IPlanck Collaboration et al.l i201lh for 
th e solar neighborhood , and comparable to the value of 0.43 found 
bv IParadis et al.l j2012ll for outer Galaxy of |h| > 10 °. Our results 
are als o in qualitative agre eme nt with the analysis bviGrenier et al.l 
( l2005ll . [Abdo et all J201oh and lAckermann et alj ^201 ih . who find 
that the DG component amounts to 40—400 per cent of the CO- 
traced mass in local small mo lecular clouds, such a s the Cepheus, 
Polaris, Chamaleon clouds etc. lMadden et al.l l ll997h find additional 
molecular mass in [C ii]-emitting regions equivalent to at least 100 
times the mass of CO traced molecular gas in irregular galaxies. 
The values of f pG in the current work are c onsistent with the the¬ 
oretical work of I Papadoj toul£s_e/aL |200^, but much larger than 
the value modelled by IWolfire et alJ ( |2010[) . /dg = 0.3, i.e. the 
DG mass is only about 43 percent of the CO traced molecular 
mass. The reason could be that the clouds we studied are all off 
the Galactic plane while IWolfire et alj ( l2010h consider the self- 
gravitating high-Av/ st ructures (Ay ^ 2 mag) tha t obey the Larson 
relations. By contrast. IPapadopoulos et alJd2002h and other studies 
are concentrate on the diffuse-type low-Av H 2 clouds that are more 
likely to be CO-dark, i.e. the type of molecular I SM th at is being 
studied in this work. Planck_CqlIalTOration et alj (l20Ilh have also 
pointe d out that IWolfire et alj ( l2010h adopted the ISolomon et al.l 
( Il987h value for the mean column d ensity of GMCs, w hich is 2-5 
times larger than the recent result of lllever et al.l ( l2009h . Amongst 
the different clouds, /dg decreases from ~ 0.7 in the Taurus ex¬ 
tensions E1-E3 (with a relatively low average extinction) to ~ 0.2 
in the Taurus and Orion clouds (with a relatively high extinction 
on average). A linear correlation between the DG column density 
Af(DG) and the visual extinction Ay is found for the different re¬ 
gions analyzed here (lower panel of Fig.|7](. The ratios of those two 
quantities show little variations amongst the different regions. We 
find that A(DG) = (1.8 ~ 2.4) x 10^'(Av — Ay)cm“^, where 
Ay is a constant, representing the minimum extinction required 
for the presence of DG. The constant Ay thus represents the ex¬ 
tinction of the H i-to-H 2 transition region. Theoretically, this value 
changes with the total H nuclei density n, metallicity z, FUV field 
Go and H 2 formation rate R f of the clouds jElmegreenll 1989lll993l : 
IPapadoDoulos et alj2002h . In the current work we find that its value 
increases with the average extinction of the region, from a value of 
0.2 mag for the North region to 1.1 mag for Orion 1. 


6 SUMMARY 

In this paper, we have constructed a 2D extinction map integrated 
out to a distance of 4 kpc for the footprint of XSTPS-GAC using the 
data presented in Paper I. The map has a high angular resolution of 
harcmin and a low noise level of 0.18 mag in Ay [~60mmag in 
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Figure 7. Upper panel: Distribution of excess column densities. White 
areas are regions where no suitable extinction or H i data are available, and 
regions where the CO emission is so intense such that Wco > 1 Kkms“*. 
The squai'es delineate the dilferent regions analyzed in the cun'ent work. 
The map has an angular resolution of 6 arcmin. Lower panel: Correlation 
between the visual extinction Ay and the DG column density. Pluses of 
different colours are binned median values for the different regions, whereas 
lines of different colours are linear fits to the binned median values. 


E{B — V)]. A co mparison with the wid ely used SFD dust map and 
the newly derived lSchlaflv et al.l < |2014|) reddening map shows that 
our map has a high fidelity. 

The integrated extinction map is combined with high- 
resolution H I and CO measurements to study the correlation be¬ 
tween the interstellar gas content and dust extinction within the 
XSTPS-GAC footprint. The main results are: 

(i) An average value of dust-to-gas ratio DGR = (4.15 + 0.01) x 
10~^^ mag cm^ and an average CO-to-H 2 conversion factor Xco = 
(1.72 + 0.01) X lO^^cm-^ (Kkms-^)”^ are found for the XSTPS- 
GAC footprint of jZ?! > 10°. The results are consistent with the 
findings of previous work. 

(ii) We find a factor of ~ 2 cloud-to-cloud variations in the 
values of DGR and Xco amongst the molecular clouds south of 
the Galactic plane within the XSTPS-GAC footprint. Variations 
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of DGR could partly be caused by the differences in dust prop¬ 
erties, such as the grain size distribution of the different clouds. 
The Xco values found for the molecular clouds are around ~ 
1. 5 X 10^° cm~^ (K km s~‘)~‘, in close agreement with the work 
of IPolket all 119881) . 


(iii) We find a range of DG mass f raction similar to those foun d 
bv IPlanck Collaboration et all 1201 ih and bv IParadis et all l2012l) . 
The mass of DG is about 23 and 124 per cent those of the atomic 
and the CO-traced molecular gas, respectively. The fraction of 
molecular mass of the component, /dg, is about 0.55. The “close” 
agreement between results from different analyses, which utilize 
different data set of differing spatial resolutions and sky cover¬ 
age probably indicate a “true universal” fraction of DG mass rel¬ 
ative to that of all gas, which is ~ 20 per cent on large Galactic 
scale. Amongst the different clouds, /dg decreases with increas¬ 
ing average extincti on. Our DG mass e stimates are much larger 
than that modelled in|WoMre_etalJj 201(]|) b ut quit e consistent with 
the theoretical work of lPapadopoulos et al.l ( l2002h . In addition, we 
find that the DG column density has a linear relationship with 
the visual extinction. The average DG-to-dust ratio N(DG)/Av is 
~ 2.0 X 10^' cm~^mag~'. The extinction of the H i -to-H 2 tran¬ 
sition layer varies between 0.22 - 1.10 mag, and increases with in¬ 
creasing average extinction of the cloud. 
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